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We supply fitting formulae enabling the normalization of slow-roll inflation models to the four-year 
COBE data. We fully include the effect of the gravitational wave modes, including the predicted 
relation of the amplitude of these modes to that of the density perturbations. We provide the 
normalization of the matter power spectrum, which can be directly used for large-scale structure 
studies. The normalization for tilted spectra is a special case. We also provide fitting functions for 
the inflationary energy scale of COBE-normalized models and discuss the validity of approximating 
the spectra by power-laws. In an Appendix, we extend our analysis to include models with a 
cosmological constant, both with and without gravitational waves. 
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I. INTRODUCTION 

The four-year data set from the Cosmic Background 
Explorer (COBE) satellite j|] is the last word we shall 
hear for some time concerning irregularities on the largest 
observable scales. One of the most important uses of the 
COBE data is in providing an accurate normalization of 
the power spectrum of density perturbations for a given 
set of theoretical assumptions; once normalized, one can 
then compare the theory to a wide range of observations 
of large-scale structure in the Universe. 

The aim of this short paper is to apply the techniques 
of Bunn and White to normalize slow-roll inflation 
models. As well as a spectrum of density perturbations, 
inflationary models predict a spectrum of gravitational 
waves, which can influence the large-angle microwave 
anisotropics seen by COBE. Further, in inflationary mod- 
els one predicts a generic link between the density per- 
turbations and gravitational waves, which must be taken 
into account in obtaining an accurate normalization. In 
this paper, we provide a self-contained account of how 
to predict these spectra from an inflationary model and 
quote fits to COBE for the normalization of the matter 
power spectrum, which can be directly used for large- 
scale structure studies. We also give fitting functions for 
the inflationary energy scale. 



II. SPECTRA FROM INFLATION 

We are of course unable to provide information for an 
arbitrary inflationary model. We shall not consider mod- 
els with more than one dynamical field, for which the 



calculations are extremely involved [gj and for which ini- 
tial conditions may be important, and we are also unable 
to consider open inflationary models ^] because in that 
case no-one has yet managed to compute the gravita- 
tional wave spectrum. We therefore restrict ourselves to 
the usual single-field inflation models, driven by a rolling 
scalar field 0, known generically as chaotic inflation 
This situation is a very general one, because most two- 
field models, for example hybrid inflation feature only 
a single dynamical scalar field, and models with extended 
gravity sectors can usually be brought into the Einstein 
form via conformal transformation [Q . 

Exact results for the spectra are not known for arbi- 
trary potentials V {(/)). They can be calculated analyti- 
cally via the slow-roll approximation. The accuracy re- 
quired depends on the way in which the results shall be 
used, and for normalizing to COBE it is valid to use the 
well known lowest-order results (see, e.g., Ref. [||), which 
give the density perturbation (i.e. scalar) spectrum As{k) 
and gravitational wave (i.e. tensor) spectrum A^^k) as 
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where prime indicates derivative with respect to cf), and 
the right-hand side is to be evaluated at the <j) value when 
the scale k equals the Hubble scale during inflation. The 
precise definition of the spectra will be clarified later. 

More accurate expressions than these do exist in the 
literature However, these are not necessary for dis- 
cussion of the COBE normalization, though they will be 
needed for discussion of the inflationary energy scale, as 
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discussed later. 

The slow-roll approximation is characterized by the 
smallness (relative to unity) of two parameters^ ||Tl|,|8t 
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Using Eq. (g) below, the spectral indices of the two spec- 
tra can be written in terms of these parameters 
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as can the ratio of the two spectra 
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Note this ratio is not independent of the tensor spectral 
index; the gravitational wave and density perturbation 
spectra are related due to their common origin in a single 
potential ^(0) 

All the above expressions apply at any scale k, with 
the spectral indices able to vary with scale. Since COBE 
covers a fairly restricted range of scales, the spectra pro- 
duced can be approximated by power-laws. Then we need 
to specify the amplitudes and spectral indices only at a 
single scale. It is best to choose this scale near the cen- 
ter of the COBE data, so we evaluate them at the scale 
fc* = TaoHo- Since, loosely speaking, the £-th microwave 
multipole samples scales around k = £aoHo/2, this cor- 
responds to the fourteenth multipole. Throughout, sub- 
script will indicate evaluation at this scale, and sub- 
script '0' indicates present value, here of the scale factor 
a and Hubble parameter H. 

For a specification of the COBE normalization to have 
a precise meaning, we need to take care in relating scales 
during inflation to present scales. The number of e- 
foldings N before the end of inflation at which k = aH 
is given by (see, e.g., Ref. ||]) 



N{k) = 62 - In 
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* Beware that two slightly different versions of these exist in 
the literature, depending on whether the fundamental quan- 
tity is taken to be the potential or the Hubble parameter . 
We use the former. 

^Only the spectra are related. The phases within a given 
realization are uncorrelated. 



Here 14 is the potential when k = aH, Vend is the po- 
tential at the end of inflation and pich is the energy den- 
sity immediately after reheating has completed, resuming 
standard big bang evolution. 

The appropriate point on the inflationary potential to 
evaluate the spectra is given by N^,. It depends on the 
energy scale of inflation, which itself depends on the nor- 
malization, fortunately only weakly. Once the COBE 
normalization is found, 14 and 14nd are determined in 
the context of a specific model, and the normalization 
can be iteratively improved if desired to take these val- 
ues into account. However, the reheat energy is much 
more uncertain, and so consequently iV* is not normally 
specified very accurately. Often, iV* is taken to be 60 or 
50. 

To locate the (f)-vahie when fc* — aH during inflation, 
one simply carries out the integral 
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where 0ond could be calculated numerically but is nor- 
mally given to adequate accuracy by the breakdown of 
the slow-roll conditions, taken as ecnd = ij^ Having lo- 
cated (/)*, one immediately gets the slow- roll parameters 
and hence the spectral indices at that scale. The tensor 
spectral index nx, which is the hardest thing to directly 
observe, can be eliminated through its relation to the ra- 
tio >1|,/Ag. To indicate the amount of tensors, we define 
a quantity r by pl|,p| 



r = 12.4 
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which measures, in the matter-dominated and Sachs- 
Wolfe approximations, the relative importance of grav- 
itational waves and density perturbations in contribut- 
ing to the relevant microwave multipole, in this case the 
fourteenth. Henceforth, n and tit will also be assumed 
to be evaluated at k*. 



III. FITTING TO THE FOUR- YEAR COBE DATA 

A. Normalization of the power spectrum 

Large-scale structure studies require the normalization 
of the present-day power spectrum. We precisely define 
our notation for the initial spectrum Aa,(k) here. In a 
critical-density universe, it is related to the rms fluctu- 
ation per logarithmic fc-interval A'^{k), and to the usual 
power spectrum P{k), both at the present epoch, by 



*In models where inflation doesn't end by steepening of the 
potential, such as hybrid inflation (pcnd is given by a dif- 
ferent condition such as an instability condition. 
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(10) 

where T{k) is the usual transfer function, normahzed to 
unity on large scales. The variance of the density field, 
smoothed on scale R, is given by 

/ ^A^{k)W'{kR), (11) 

where the smoothing function W{kR) tends to unity at 
small k. The observables related to are discussed in 
Refs. @. 

The fitting to COBE is described in Ref. |]. Rather 
than use Eq. (^ to set the relative normalization of 
scalars and tensors, we use the more accurate expression 

This relation comes from a full higher-order calculation, 
so higher-order expressions can also be used if desired to 
compute n and nx from a given inflation model, along 
the lines of Refs. |^,|l^. Using this, r is related to the 
spectral indices by 

r = -6.2nT(l-^ + (7i-l)) . (13) 

We specify the normalization of the density pertur- 
bations, following Ref. at the present Hubble scale 
k = QqHq^ and define 

6H = AsiaaHo). (14) 

By focusing on the normalization at such a large scale, 
we obtain a result which is independent, to excellent ac- 
curacy, of the choice of cosmological parameters, such as 
the present Hubble constant and the nature of the dark 
matter. The one exception is a nonzero cosmological 
constant, which does affect the large-angle anisotropics; 
we generalize our results to that case in the Appendix. 

Our main result is a fitting function for the COBE 
normalization, which is accurately represented by 

fe(n,.)^ 1.91x10-^ '''!!!?.r^^- (15) 

Vl + O.Yor 

The Icr observational error is 7%. Within the region 
0.7 < n < 1.3 and —0.3 < ut < (the latter corre- 
sponding to < r < 2), the fit is good to within 1.5% 



everywhere. The change from varying other cosmological 
parameters (except the density parameter) is within 4% 
for reasonable variations. In addition, there is a system- 
atic uncertainty of ~ 3% associated with the process by 
which the COBE maps are made. Combining all of these 
uncertainties in quadrature, we believe that a realistic 
estimate of the uncertainty in (5h is 9% at la. 

The terms in Eq. ( [l5| ) have a simple interpretation. 
The numerical prefactor is the result for a scale-invariant 
density perturbation spectrum. The n term represents 
the pivot point of the COBE data; it guarantees that 
COBE normalized spectra at fixed r cross at fcpivot = 
e^-^^aoHo. This number actually corresponds to about 
the fifteenth multipole (and in fact a purely scalar fit 
even prefers the sixteenth), but the tensors give greater 
weight to the lower multipoles, making £ = 14 a better 
overall choice for the pivot. The pivot point is at higher 
£ in the four-year data than in the two-year data, since 
the low multipoles were already cosmic-variance limited 
and the higher i have improved signal-to-noise ratio. 

The final term in Eq. ( |l5|) is the reduction in amplitude 
due to the tensors. It is interesting that the tensor term 
has a coefficient of only 0.75, since the definition of r was 
intended to make that coefhcient close to unity. How- 
ever, the factor 12.4 in Eq. (||) was computed ||ll| using 
the fully matter-dominated Sachs- Wolfe approximation 
for both the tensor and scalar spectra. The dominant 
correction to this is the effect on the tensor spectrum 
from the universe not being perfectly matter dominated 
at last scattering ||l^ (about twenty percent), with the 
start of the rise to the acoustic peak from the density 
perturbations contributing another five percent. These 
corrections have been noted in papers concerned with 
'cosmic confusion' ||l^. That the coefRcient is only 0.75 
means that papers using the original analytic argument 
have somewhat over-estimated the amount by which ten- 
sors reduce the power spectrum normalization. 

B. The inflationary energy scale 

For a given inflationary model, the COBE normaliza- 
tion fixes the energy scale of inflation at that time, V^,. 
The fitting function can therefore be inverted to supply 
the inflationary energy scale [ p"5|JT^ . One would like to 
use Eq. (|l|) in order to do this, but in fact this equation is 
not always very accurate at determining the amplitude, 
since the slow-roll approximation is not necessarily as ac- 
curate as the observations and fits we have discussed. We 
therefore use the next-order version of this equation, as 
derived by Stewart and Lyth Q , which is^ 



^It is irrelevant that this differs from the scale at which we 
evaluated the spectra to do the normalization. The two scales 
are simply related, As{aoHo) = 7^""^i(7ao//o). 

**For the record, we take pure cold dark matter with h — 

0.75 and nbaryon/l^ = 0.0125. 



^'To derive this from their results, some further relations are 
necessary which can be found in Ref. |lO[|. 
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It is vital to note that we did not have to use this for the 
normalization, because the amplitude correction in the 
prefactor is typically almost constant across the COBE 
scales (that is, e and r] hardly vary). It therefore cancels 
out when one does the COBE normalization for large- 
scale structure. 

We specify the energy at fc* , the place where the slow- 
roll parameters were evaluated. Since A^(fc*) = 7" ^^(5^, 
substituting in the fitting function and carrying out a 
small parameter expansion gives 

V,^^^ = (5.4 X 10-3 mpi) tj^ (1 - 3.2e, -I- 0.5?7,) , (17) 
= (6.6 X 10^^ GeV) tj^ (1 - 3.2e, + 0.5?7,) . 

Even if one does not have a specific inflation model in 
mind, it is possible to use this to obtain an upper bound 
on the energy density at the end of inflation by imposing 
some assumptions about how inflation will end 



C. On the validity of the power-law approximation 

It is possible to test the validity of approximating the 
spectra by power-laws. We'll concentrate on the density 
perturbations. To a first approximation, variation in the 
spectral indices is driven by the difference between n—\ 
and nx, but one should consider the full slow-roll formula 
for that variation at fc* which is llTtl 
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The COBE data extend only for slightly more than one 
log interval in either direction about the central point, so 
unless this number is greater than, say, a few hundredths, 
the power-law approximation for COBE will be excellent. 
Kosowsky and Turner |l^ evaluated this for a range of 
inflation models without flnding a value anywhere near 
this large. 

It is possible that the spectra may be well approx- 
imated by power-laws on the COBE scales, but not 
across the much wider range corresponding to future mi- 
crowave anisotropy observations and to large-scale struc- 
ture [|l^,0 . In that case, one simply uses our results with 
the appropriate approximate power-law at the COBE 
scales. 



IV. A WORKED EXAMPLE: THE QUADRATIC 
POTENTIAL 

The simplest inflationary model is chaotic inflation |^] 
with a quadratic potential 1^(0) = ra^^P' /2. For this 
potential, the slow-roll parameters are 



1 TOpi 

47r dp- 



and hence 



(19) 



TOpi/\/47r- For definiteness, we take 



iV* = 60, and from Eq. (|^) we flnd 
this model one predicts 



3.10 mpi. So for 



n = 0.967 



nx = -0.017 



r = 0.10 



(20) 



The normalization of the matter power spectrum, from 
Eq. (^, is therefore = 1.91 x 10"^. Ahhough this 
model is very close to scale-invariant limit, the difference 
is still non-negligible. For example, if one computes the 
variance at d>h~^ Mpc, denoted erg, it is reduced, relative 
to scale-invariance with no gravitational waves, by 10%, 
where 6% is due to the tilt and 4% due to the gravi- 
tational waves. Many other inflation models give much 
larger corrections than this. 

The inflationary energy scale corresponding to this, 
from Eq. (^, is 



y^'' 1.6 X IQ-^ mpi = 2.0 X 10^^ GeV , 



(21) 



corresponding to m — 1.2 x lO^^mpi. 

Finally, the scale-dependence of the spectral index is 
dn/d\nk\f: = —5 x 10^'*, which is completely negligible 
for COBE. 



V. SUMMARY 

One of the most important uses of the COBE data is 
to normalize the matter power spectrum used in large- 
scale structure studies. We have shown that it is pos- 
sible to condense this information into a single fltting 
function, covering not just the case of tilted perturba- 
tion spectra but also including the spectrum of gravita- 
tional waves that inflation predicts. A further extension 
to models with a cosmological constant is given in the 
Appendix. Within the slow-roll paradigm, the amount 
of gravitational waves, parameterized by r, is completely 
independent of the density perturbation slope n, though 
the scale-dependence of the gravitational wave spectrum 
is then predicted ||ll|,|| . The fitting functions can also be 
used for the case without gravitational waves, simply by 
setting r equal to zero. 
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APPENDIX A: GENERALIZATION TO MODELS 
WITH A COSMOLOGICAL CONSTANT 

The dynamics of inflation are insensitive to whether or 
not there is a present cosmological constant A, provided 
the spatial geometry is kept flat. We now provide the 
generalization of our expressions to the case where A 0. 

We need a generalization of Eq. to account for the 
change in the growth of perturbations in a low-density 
universe, and the relation between the curvature pertur- 
bation and the matter power spectrum. If Su continues to 
indicate the present power spectrum, and the initial 
perturbation spectrum as before, then 



6h 



g(»o) 

fin 



AsiaoHo) , 



(Al) 



and the right hand side of Eq. is multiplied by 

f/^(rio)/f^o- Here g{fl) is the growth suppression factor, 
which is accurately fit by |l9| 



m = 2^ 



1 

70 
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(A2) 



A generalization of the factor 0.75 multiplying the 
gravitational wave term in the fitting function for Su is 
also needed. By directly evaluating the radiation power 
spectra for scale-invariant initial spectra with As = At = 
1, we find the ratio of contributions to C14 is fit to within 
a percent by 



^ ^ /(f^o) ^ 0.75 

'-'14 



0.13 r^i. 



(A3) 



with Q\ = 1 — fio- 

First considering the case with no tensors, a good fit 
is obtained by using an Qq dependence plus a cross-term 
between n = n — 1 and The formula 

Sn{n, no) = 1.91 x 10"^ exp [-1.01 H] o.8o-o.o5i„Oo 
X [l + 0.18Hf7A] , (A4) 

holds within 2.5% for n as before and 0.2 < < 1. 

We can extend this to include tensors using Eq. (A3) 
and the addition of an extra cross-term which vanishs 
in the critical-density case and the tensorless case. This 
gives 

6u = 1.91 X 10-^ e xp [-1.01 n] o.80-o.05inf2o 

X [1 + 0.18 nf^A - 0.03 rf^A] • (A5) 

For the parameter ranges 0.7 < n < 1.3, —0.3 < nx < 
and 0.2 < flo < 1, this fit is within 1% almost every- 
where, and always within 2.5%. The relation between 
riT and r is still given by Eq. ([l^). 
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